A 147  011 


AFOSR-TR-  84  -0  87  8 


SC5361.6SAR 
Copy  No.*^  3 


NONLINEAR  WAVE  PROPAGATION  STUDY 


SEMI-ANNUAL  TECHNICAL  REPORT  NO.  2  FOR  THE  PERIOD 
December  1,  1983  through  May  31,  1984 


DTIC 

ELECTEI 


NOV  1 


CONTRACT  NO.  F49620-83-C-0065 
DARPA  ORDER  NO.  4400 
PROGRAM  CODE:  3A10 


Prepared  for 

Air  Force  Office  of  Scientific  Research 
Building  410 
Bolling  AFB,  DC  20332 


Principal  Investigator 
B.R.  Tittmann 
(805)  498-4545 


Sponsored  by  ' 

Defense  Advanced  Research  Projects  Agency  (DoD) 
DARPA  Order  No.  4400 

Monitored  by  NP  Under  Contract  NO.  F49620-83-C-0065 


The  views  and  conclusions  contained  in  this  document  are  those  of  the  authors 
and  should  not  be  interpreted  as  necessarily  representing  the  official  policies, 
either  expressed  or  implied,  of  the  Defense  Advanced  Research  Projects  Agency 
or  the  U.S.  Government. 


approved 


Rockwell  International 

Science  Center 


sy  /o  ay  on 


1»  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED  _ 


2*  security  classification  authority 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


^2b  OCCLASSIPICATION/OOVMNGRAOING  SCHEDULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

SC5361 . 6SAR 

64  NAME  of  PERFORMING  ORGANIZATION 

Sb.  OFFICE  SYMBOL 

Rockwell  International 

Science  Center 

Of  applicable) 

6c  AOORESS  (City.  State  and  ZIP  Code) 

1049  Camino  Dos  Rios 

PO  Box  1085 

Thousand  Oaks,  CA  91360 

to  NAME  OF  FUNDING/SPONSORING 

8b.  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable ) 

AFOSR/NP 

NP 

8c  aODRESS  (City.  State  and  ZIP  Code) 

Btiidling  410 

Foiling  AFB  DC  20332-6448 

3.  DISTRIBUTION/AVAILABILITY  OF  REPORT- 

Approved  for  public  release; 
Distribution  unlimited 


(^ORGANIZATION  REPORT  NUMBERIS) 

SR-TR-  8  4.0  878 


AFOSR/NP 


7b.  ADDRESS  (City.  State  and  ZIP  Code ) 

Building  410 

Bolling  AFB  DC  20332-6448 


F49620-83-C-0065 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 
ELEMENT  NO. 


I’oNLtirffi' w^e'^roTa'gatton  study 


62714E 


PROJECT 

NO. 

t  O 


WORK  UNIT 
NO 


12,  PERSONAL  AUTMOR(S) 

B.  R.  Tittman 


13*.  TYPE  OF  REPORT 

SEMI-ANNUAL  TECH 


13b.  TIME  COVERED 


14.  DATE  OF  REPORT  (Yr..  Mo..  Day)  I  15.  PAGE  COUNT 


e!S££?9sfi3 


31  May  84 


COSAT  I  COOES 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


tf  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number* 

In  this  report  the  results  of  a  study  of  the  response  of  Westerly  granite 
to  sinusoidal  loading, are  presented.  The  amplitude  of  transition  from 
linearity  to  non-linearity  can  be  defined.  Results  are  compared  to  previous 
studies. 


20  0»STRiBUTION/A  V  At  LABILITY  Of  ABSTRACT 
^UNCLASSiFiEO/UNLIMiTEO^B  SAME  AS  RPT.  Q  DTlC  USERS  D 


23t.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Henry  R.  Radoski 


..$0.FORMH73,83APR 


T  21.  ABSTRACT  SECURITY  CLASSIFICATION 

3  DT .  c  USE  RS  □  UNCLAS  S I  FI  ED 


EDITION  OF  1  JAN  73  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


22b  TELEPHONE  NUMBER  22c  OFFICE  SYMBOL 

(include  Area  Code} 

202/767-4906  NP 


Rockwell  International 

Science  Center 


SC5361.6SAR 


TABLE  OF  CONTENTS 


Page 


1.0  INTRODUCTION .  1 

2.0  EXPERIMENTAL  STUDIES .  2 

3.0  DISCUSSION .  13 

4.0  SUMMARY  AND  CONCLUSIONS .  14 

5.0  REFERENCES .  15 


A1H  TORCS  OJ7IC5  W  SCTESTITTC  r  STS?’ 

NOTICE  0?  rVJtW-.TTiT:  TO  I>?tC 

This  ••  ''  *  1S 

approve  .;  :  or  *v.  .  ’  •'  1  :;f  *“12* 

Distributio.-  V.’  lilted. 

MirrHEW  J.  KSSF'AR 

Chief.  leehalcal  Inforaatloa  Dlrlelon 


i  i  1 

C6209A/md 


•V -TN 


r V/.  ’A"1?  V  «.'*/  i,'  *.111 1.<  l  '  ■  •  q 


Rockwell  International 

Science  Center 


SC5361.6SAR 


LIST  OF  FIGURES 


Figure  Page 

1.  Definition  of  strain  amplitude  in  flexure  and  torsion .  3 

2.  Attenuation  and  log  RMS  driving  transducer  voltage 


plotted  as  a  function  of  strain  amplitude  for  torsional 
vibrations  in  dry  Westerly  granite  at  0  MPa  effective 


pressure . 3 

3.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
torsional  vibrations  in  dry  Westerly  granite  at 

1.7  MPa  effective  pressure .  4 

4.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
torsional  vibrations  in  dry  Westerly  granite  at 

3.4  MPa  effective  pressure .  4 

5.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
torsional  vibrations  in  dry  Westerly  granite  at 

6.8  MPa  effective  pressure .  5 

6.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
torsional  vibrations  in  dry  Westerly  granite  at 

34  MPa  effecitve  pressure .  5 

7.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
flexural  vibrations  in  dry  Westerly  granite  at 

0  MPa  effective  pressure .  6 

8.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
flexural  vibrations  in  dry  Westerly  granite  at 

1.7  MPa  effective  pressure .  6 


K.v 


9. 


Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for  flexural 
vibrations  in  dry  Westerly  granite  at  3.4  MPa  effec¬ 
tive  pressure . 


v  1  v 

;•  C6209A/md 


r«*w‘.  «r. •*.«  ,  •* 


Rockwell  International 

Science  Center 


B 


tv 


; 


¥ 


SC5361.6SAR 


10.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
flexural  vibrations  in  dry  Westerly  granite  at 

6.8  MPa  effective  pressure .  7 

11.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
flexural  vibrations  in  dry  Westerly  granite  at 

34  MPa  effective  pressure.... .  8 

12.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
torsional  vibrations  in  wet  Westerly  granite  at 

0  MPa  effective  pressure .  8 

13.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
torsional  vibrations  in  wet  Westerly  granite  at 

1.7  MPa  effective  pressure .  9 

14.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
torsional  vibrations  in  wet  Westerly  granite  at 

3.4  MPa  effective  pressure .  9 

15.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
torsional  vibrations  in  wet  Westerly  granite 

6.8  MPa  effective  pressure .  10 

16.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
torsional  vibrations  in  wet  Westerly  granite  at 

34  MPa  effective  pressure .  10 

17.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
flexural  vibrations  in  wet  Westerly  granite  at 

0  MPa  effective  pressure .  11 

18.  Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
flexural  vibrations  in  wet  Westerly  granite  at 

1.7  MPa  effective  pressure .  11 


v 

C6209A/md 


Rockwell  International 

Science  Center 


SC5361.6SAR 

Attenuation  and  log  RMS  driving  transducer  voltage 


plotted  as  a  function  of  strain  amplitude  for 

flexural  vibrations  in  wet  Westerly  granite  at 

3.4  MPa  effective  pressure .  12 

Attenuation  and  log  RMS  driving  transducer  voltage 

plotted  as  a  function  of  strain  amplitude  for 

flexural  vibrations  in  wet  Westerly  granite  at 

6.8  MPa  effective  pressure . . .  12 

Attenuation  and  log  RMS  driving  transducer  voltage 
plotted  as  a  function  of  strain  amplitude  for 
flexural  vibrations  in  wet  Westerly  granite  at 

34  MPa  effective  pressure .  13 


vi 

C6209A/md 


.V-V.VA  .  .  V.V  \ 


Ml 


KtfCi 


•.-  %• «.  ■  . 


9 


Rockwell  International 

Science  Center 


SC5361.6SAR 


1.0  INTRODUCTION 


At  the  present  time  there  still  exists  an  element  of  uncertainty 
regarding  the  definition  of  an  effective  radius,  which  can  be  used  to  separate 
the  nonlinear  near-field  regime  adjacent  to  an  explosion  from  the  linear  far- 
field  seismic  regime.  Using  data  sets  from  the  COWBOY  experiments  and  the 
SALMON  test  Trulio  (1978,  1981)  studied  the  decay  of  peak  particle  velocity  as 
a  function  of  scaled  distance  in  dome  salt.  He  observed  that  decay  rates  are 
significantly  faster  than  can  be  explained  assuming  perfectly  elastic  behav¬ 
ior,  and  concluded  that  material  behavior  was  not  perfectly  elastic,  and  prob¬ 
ably  nonlinear,  even  at  strains  as  low  as  those  encountered  at  the  most  dis¬ 
tant  instruments.  Larson  (1982)  examined  near-field  particle  motions  in 
pressed  salt  in  the  laboratory  using  small-scale  chemical  explosions.  He 
demonstrated  near-linear  superposition  of  wave  forms  even  at  strains  in  excess 
of  10"4.  Yet,  the  Q  values  calculated  by  Larson  (12.5  near  10-3  strain  and 
24.9  near  6  10-4  strain)  depended  on  amplitude,  indicating  nonlinear 

behavior.  Furthermore,  these  values  are  low  compared  with  those  measured  by 
Tittmann  (1983a)  using  the  same  material  except  at  much  lower  amplitudes.  The 
work  of  Tittmann  (1983a,  1983b)  demonstrates  that  nonlinear  effects,  as  re¬ 
flected  in  an  amplitude-dependent  Q,  can  persist  to  very  low  levels  of  strain, 
even  below  10-6  in  some  cases.  More  recently,  Minster  (1982)  and  Minster  and 
Day  (1984)  re-examined  the  COWBOY  data  and  concluded  that  the  observed  decay 
rates  can  only  be  explained  using  an  amplitude  dependent  Q.  However,  Burdick 
et  al  (1984a)  have  argued  that  it  is  possible  to  model  a  seismic  source  func¬ 
tion  for  the  Amchitka  tests  assuming  linear  behavior  in  the  near  field  just 
outside  the  spall  zone  (approx.  700-1200  m/kt^3).  Furthermore,  Burdick  et  al 
(1984b)  argued  that  the  same  model  can  be  used  to  predict  the  first  vertical 
pulse  arrival  and  rise  times  even  within  the  spall  zone.  They  used  the  con¬ 
cept  of  a  compressional  elastic  radius,  which  may  in  fact  be  considerably  less 
than  a  tensional  elastic  radius,  which  must  extend  at  least  as  far  as  the 
outer  limits  of  the  spall  zone. 
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In  this  report  we  present  the  results  of  a  recent  study  of  the 
response  of  Westerly  granite  to  sinusoidal  loading.  The  amplitude  of  tran¬ 
sition  from  linearity  to  nonlinearity  can  be  defined,  and  results  are  compared 
to  previous  studies. 


2.0  EXPERIMENTAL  STUDIES 

The  techniques  used  to  obtain  the  data  presented  in  this  report  are 
the  same  as  those  described  in  an  earlier  report  (Tittmann,  1983a).  Briefly, 
a  cylindrical  specimen  of  Westerly  granite,  13  cm  in  length  and  0.71  cm  in 
radius,  was  clamped  into  an  apparatus  which  enabled  it  to  be  vibrated  in  both 
flexural  and  torsional  modes.  In  this  particular  study  the  resonant  frequency 
was  between  500  and  550  Hz  for  all  measurements.  The  Q  of  the  bar  could  be 
calculated  from  the  ratio  of  the  resonant  frequency  to  the  bandwidth  of  the 
resonance  curve.  The  voltage  of  the  electromagnetic  transducer,  used  to  drive 
the  resonance,  could  be  adjusted  so  as  to  vary  the  amplitude  of  vibration 
(vibration  amplitudes  for  torsional  and  flexural  vibrations  are  defined  in 
Fig.  1).  Thus,  attenuation  (Q"1)  and  RMS  driving  voltage  can  be  measured  as  a 
function  of  vibration  amplitude  in  the  sample.  Provided  that  the  sample 
behaves  linearly,  the  attenuation  should  be  independent  of  vibration  amp¬ 
litude,  and  the  RMS  driving  voltage  should  be  proportional  to  vibration  amp¬ 
litude  (c.f.  Tittmann,  1983a).  Measurements  were  made  in  both  wet  and  dry 
Westerly  granite  at  effective  pressures  from  0  to  34  MPa. 

The  results  of  this  study  are  shown  in  Figs.  2-21.  The  results  are 
qualitatively  consistent  with  the  results  of  previous  studies  (Tittmann, 

1983a)  in  that  linear  behavior  is  observed  at  all  low  amplitudes,  with  a 
transition  for  nonlinear  behavior  at  high  strain.  The  transition  amplitude  is 
near  10’6  in  all  cases.  Including  both  torsion  and  flexure,  increasing  very 
slightly  with  increasing  effective  pressure. 
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Fig.  1  Definition  of  strain  amplitude  in  fluxure  and  torsion. 
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Fig.  2  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  torsional  vibrations 
in  dry  Westerly  granite  at  0  MPa  effective  pressure. 
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Fig.  4  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  torsional  vibrations 
In  dry  Westerly  granite  at  3.4  MPa  effective  pressure. 
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LOG  MAXIMUM  STRAIN  AMPLITUDE 

Fig.  5  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  torsional  vibrations 
in  dry  Westerly  granite  at  6.8  MPa  effective  pressure. 
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Fig.  6  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  torsional  vibrations 
in  dry  Westerly  granite  at  34  MPa  effective  pressure. 
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LOG  MAXIMUM  STRAIN  AMPLITUDE 

Fig.  7  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
in  dry  Westerly  granite  at  0  MPa  effective  pressure. 
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Fig.  8  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
in  dry  Westerly  granite  at  1.7  MPa  effective  pressure. 
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LOG  MAXIMUM  STRAIN  AMPLITUDE 

Fig.  9  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
in  dry  Westerly  granite  at  3.4  MPa  effective  pressure. 
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Fig.  10  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
in  dry  Westerly  granite  at  6.8  MPa  effective  pressure. 
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Fig.  11  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
in  dry  Westerly  granite  at  34  MPa  effective  pressure. 
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Fig.  12  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  torsional  vibrations 
in  wet  Westerly  granite  at  0  MPa  effective  pressure. 
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LOG  MAXIMUM  STRAIN  AMPLITUDE 

Fig.  13  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  torsional  vibrations 
in  wet  Westerly  granite  at  1.7  MPa  effective  pressure. 
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Fig.  14  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  torsional  vibrations 
in  wet  Westerly  granite  at  3.4  MPa  effective  pressure. 
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LOG  MAXIMUM  STRAIN  AMPLITUDE 

Fig.  15  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  torsional  vibrations 
in  wet  Westerly  granite  at  6.8  MPa  effective  pressure. 
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Fig.  16  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  torsional  vibrations 
In  wet  Westerly  granite  at  34  MPa  effective  pressure. 
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Fig.  17  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
in  wet  Westerly  granite  at  0  MPa  effective  pressure. 
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Fig.  18  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
In  wet  Westerly  granite  at  1.7  MPa  effective  pressure. 
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LOG  MAXIMUM  STRAIN  AMPLITUDE 


Fig.  19  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
in  wet  Westerly  granite  at  3.4  MPa  effective  pressure. 
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Fig.  20  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
in  wet  Westerly  granite  at  6.8  MPa  effective  pressure. 
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LOG  MAXIMUM  STRAIN  AMPLITUDE 

Fig.  21  Attenuation  and  log  RMS  driving  transducer  voltage  plotted 
as  a  function  of  strain  amplitude  for  flexural  vibrations 
in  wet  Westerly  granite  at  34  MPa  effective  pressure. 

3.0  DISCUSSION 

In  this  phase  of  the  study,  as  in  those  presented  earlier,  the  test 
specimen  was  subjected  to  a  long  train  of  sinusoidal  vibrations.  For  flexural 
vibrations  the  test  material  spent  half  of  its  cycle  primarily  in  tension  and 
half  of  its  cycle  primarily  in  compression.  Therefore,  using  this  type  of  mea¬ 
suring  technique  with  harmonic  vibrations  it  is  not  possible  to  establish 
whether  the  nonlinearity  observed  in  the  flexural  mode  data  is  due  to  re¬ 
laxation  in  compression,  relaxation  in  tension,  or  both.  Recognizing  the  fact 
that  rocks  are  significantly  stronger  in  compression  than  In  tension  it  is 
reasonable  to  expect  that  the  nonlinearity  is  primarily  a  tensional  phenome¬ 
non,  but  this  remains  to  be  established  experimentally.  More  data  on  the  con¬ 
stitutive  properties  of  rocks  at  modest  amplitudes  (10“^  to  10“®)  are  re¬ 
quired.  An  Ideal  study  would  Include  measurements  In  both  uniaxial  compres¬ 
sion  and  uniaxial  tension.  Other  issues  regarding  the  applicability  of  the 
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resonating  bar  technique  for  simulating  non-linear  pulse  loading  are  discussed 
in  the  previous  report  (Tittman,  1983b),  and  will  not  be  repeated  here. 

It  is  quite  reasonable  to  assume  that  the  amplitude  of  transition 
from  linear  to  nonlinear  behavior  detected  with  the  resonating  bar  apparatus 
represent  a  lower  limit  on  the  corresponding  transition  amplitude  for  a 
seismic  pulse  propagating  in  the  near-field  of  an  explosion.  The  accumulated 
damage  caused  by  subjecting  the  specimen  to  repeated  cycles  of  high  amplitude 
vibration  (Tittmann,  1983b)  would  have  the  effect  of  making  the  transition 
amplitude  appear  relatively  low  in  laboratory  experiments  compared  with  the 
field.  Furthermore,  the  material  near  the  explosion  is  subjected  to  lower 
levels  of  tensional  stress  than  compressional  stress,  and  tensional  non¬ 
linearity  may  occur  at  lower  stress  levels  than  compressional  nonlinearity. 

Measurements  of  the  constitutive  properties  of  materials  in  both 
tension  and  compression  are  critical  and  will  be  included  in  future  studies. 


4.0  SUMMARY  AND  CONCLUSIONS 

'  ^  :  ' 

Nonlinear  effects^  have  been  observed  in  a  cylindrical  test  specimen 
of  Westerly  granite  which  was  subjected  to  both  flexural  and  torsional  modes 
of  resonant  vibration.  Nonlinear  effects  in  shear  are  observed  when  the 
shearing  strain  exceeds  approximately  10"®,  increasing  slightly  with 
increasing  effective  pressure.  Nonlinear  effects  in  flexure  are  also  observed 
when  the  extension/compression  strain  exceeds  10”®,  also  increasing  very 
slightly  with  increasing  effective  pressure.  These  transition  amplitudes 
probably  represent  a  lower  limit  on  the  amplitude  of  transition  from  linear  to 
nonlinear  behavior  for  the  primary  elastic  pulse  propagating  in  the  near-field 
of  an  explosion.  Using  these  measurements  it  is  not  possible  to  determine 
whether  the  nonlinear  effects  observed  In  flexure  are  primarily  an  extensional 
or  compressional  feature.  This  is  a  critical  issue  and  in  future  studies  we 
will  attempt  to  examine  the  constitutive  properties  of  rocks  in  both  tension 
and  compression  under  moderate  to  low  amplitude  conditions. 
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